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We experimentally investigate the dielectric response of the low-dimensional gapped quantum 
magnet Cu2Cl4-HgC4S02 near a magnetic field-induced quantum critical point, which separates the 
quantum-disordered and helimagnetic ground states. The observed magnetocapacitive effect origi¬ 
nates from an improper ferroelectric nature of the transition, which itself is perhaps one of the best 
known realizations of Bose-Einstein condensation of magnons. Despite that, we find that the mag¬ 
netocapacitive effect associated with the transition exhibits huge and very unusual anharmonicities. 
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The phenomenon of magnetic Bose-Einstein conden¬ 
sation (BEC) is one of the most fascinating analogies 
between the properties of quantum magnets and bosonic 
systems mm- It corresponds to spontaneous long-range 
ordering induced in gapped quantum paramagnets by 
external magnetic fields. Since conventional BEC is a 
spontaneous breaking of 1 /( 1 ) gauge symmetry, a cru¬ 
cial requirement for the magnetic analog is the pres¬ 
ence of isomorphous SO{2) axial spin rotation symme¬ 
try in the disordered phase. Unfortunately, symmetry¬ 
breaking anisotropic interactions that are always present 
in real crystalline materials preclude the applicability of 
the BEC concept close to the quantum critical point. The 
only exception would be the case of incommensurate he¬ 
limagnetic order. Here, the required SO{2) symmetry is 
robust, as it represents a sliding of the magnetic spiral 
with respect to the crystal lattice, made “frictionless” 
by its incommensurability. There is, however, a poten¬ 
tial complication. Helimagnets are likely to be multi- 
ferroic [3]. In a spiral structure, the lack of inversion 
symmetry between the spins gives rise to electric polar¬ 
ization [1H6]. The intrinsic coupling between magnetic 
and electric degrees of freedom is a complication that 
can potentially modify the very nature of the quantum 
critical point (QCP). A few recent studies looked at rare 
cases in which multiferroicity occurs at quantum phase 
transitions (QPTs) but not much work was done 

to study the criticality of such transitions. Moreover, the 
very important case of QCPs connecting magnetically 
ordered and quantum-disordered phases remains largely 
unexplored in the context of multiferroic physics. The 
two questions to be addressed are as follows: (i) Can this 
type of transition be described as BEC and (ii) how do 
the dielectric properties behave at and beyond the QCP? 

On the experimental side, some progress has occurred 
very recently with the discovery of the S' = 1/2 quasi-one- 
dimensional quantum antiferromagnet Cu 2 Cl 4 -H 8 C 4 S 02 
(also known as SUI-CU 2 CI 4 ) This material fea¬ 

tures a magnetic-field-induced QPT from a quantum- 
disordered to a helimagnetic state. However, previous 
studies have portrayed this transition as being more 
complex than a simple magnetic BEC. Dielectric spec¬ 
troscopy experiments [14] pointed to the possibility of a 
dielectric susceptibility divergence at T = 0. At the same 


time, the shape of the phase boundary and the order pa¬ 
rameter exponent were found to be inconsistent with the 
BEC paradigm [iQiiig. This indicated the possibility of 
an unconventional order parameter involving both spin 
and charge degrees of freedom. In the present Rapid 
Communication, we show that this is not the case, and 
that electric polarization in SUI-CU 2 CI 4 is not a true or¬ 
der parameter of the transition. We show that the QCP 
is a “protected” BEC of magnons, and is arguably the 
only clean realization of such among all known materi¬ 
als. While electric polarization plays a dependent role, 
its behavior is far from trivial. Near the magnetic QCP 
in SUI-CU 2 CI 4 , we find spectacular anomalies in the non¬ 
linear dielectric response. 

SUI-CU 2 CI 4 belongs to a family of insulating metaloor- 
ganic Heisenberg spin systems. Its magnetic properties 
are due to S' = 1/2 Cu^+ ions. The spins are an- 
tiferromagnetically coupled into one-dimensional struc¬ 
tures, which can be described as highly frustrated four- 
leg spin tubes, running along the c axis of the triclinic 
structure (see Refs. mm for more details). This lay¬ 
out, with an even number of “legs”, is responsible for 
the nonmagnetic quantum-disordered ground state, and 
a gap A ~ 0.52 meV in the spin excitation spectrum. 
Due to a geometric frustration of the exchange interac¬ 
tions, the dispersion minimum for the triplet of lowest- 
energy S' = 1 excitations occurs at an incommensurate 
wave vector. As a result, in an external magnetic field 
He = A/gfiB ~ 3.7 T for H||b, SUI-CU 2 CI 4 undergoes 
a soft-mode ordering transition with a propagation vec¬ 
tor Q = (—0.22,0,0.48) [T2l[T3]. The structure of the 
ordered high-held phase has been resolved by neutron 
diffraction m- The spin components transverse to H 
form a spiral such that (Sj_(r)) = Si cos(Qr)-|-S 2 sin(Qr) 
with Si T S 2 . This spiral arrangement lacks an inver¬ 
sion symmetry and hence induces an electric polarization 
P oc [[Si X S 2 ] X Q] mis]. To date, this polarization has 
not been directly measured, and is likely to be extremely 
small. Instead, previous studies of Schrettle et al. [14] 
probed the corresponding dielectric permittivity e. The 
apparent divergence of this quantity at the phase transi¬ 
tion was taken as a sign of the critical huctuations of P, 
which would imply that the transition is different from 
simple ordering of localized spins. 
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FIG. 1. (Color online) Magnetocapacitive effect in Sul- 
CU2CI4. The anomalous contribution to the sample capaci¬ 
tance AG is shown as a function of temperature and magnetic 
field, (a) Field scans at constant temperature. The peak po¬ 
sitions and amplitudes are projected onto the AG — H and 
AC — T planes (yellow diamonds). The highlighted H — T re¬ 
gion in close vicinity to the QCP is the domain of temperature 
scans, (b) Temperature scans. The step in field is 0.002 T. 
The extrapolated value of the zero-temperature critical field 
is marked, (c) Measured dependence of the amplitude of the 
dielectric anomaly on the probing voltage at T = 500 mK. 


The main limitation of that study was in that it only 
probed temperatures above 1 K, never closely approach¬ 
ing the QCP at T = 0. In the present work we have over¬ 
come this technical difficulty by combining a capacitance 
bridge setup similar to that used in Ref. m with a ^He- 
“^He dilution refrigerator [15]. What is in fact measured 
in our experiments is the capacitance of a plate capacitor 
with a SUI-CU 2 CI 4 sample in between the plates. At low 
temperatures in the absence of magnetic field it is domi¬ 
nated by a constant background Cq ~ 3.5 pF due to the 
“normal” dielectric permittivity of SUI-CU 2 CI 4 , £ ~ 3. 
The quantity of interest is the additional capacitance 
AC{H, T), which is both temperature and magnetic field 
dependent near the QCP. It is plotted in Figs. and 
0 b) and represents the observed magnetocapacitive ef¬ 
fect, with a maximal detected Ae/e ~ 3%. 

The prominent capacitance peak seen at high tempera- 



FIG. 2. (Color online) A few representative field scans of a 
sample capacitance (top), compared to the specific heat data 
(bottom). For both measurements, H || a. 
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FIG. 3. (Color online) Phase diagram of SUI-CU 2 CI 4 for H | 
a. Small circles are the transition points determined from 
the magnetocapacitive effect and diamonds are the transition 
points found from the specific heat A-anomaly. The solid line 
is a guide for the eye. The dashed line is a fit to Eq. 0 in 
the range T < 300 mK. The inset shows a logarithmic plot of 
the measured phase boundary. 


tures is what Ref. m must have taken for a divergent e. 
In fact, at low temperatures, as the QCP is approached, 
the anomaly weakens, and transforms into a small and 
rounded step at T —3 0. The peak amplitude of AC 
decreases with the temperature almost linearly, rather 
than diverging at the QCP. A comparison between the 
capacitance peaks and A-anomalies in specific heat for 
H II a (Fig. [|) shows that the transition point corre¬ 
sponds to the inflection point of the AC{H) steep slope. 
Note that critical fluctuations are prominent in specihc 
heat well below the transition held, while the dielectric 
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contribution immediately disappears in the disordered 
phase. Moreover, as shown in Fig. [^b), there is no 
detectable change in the dielectric permittivity outside 
of the ordered phase even along the eritical trajectory 
H = He = 3.43 ± 0.01 T [T^. We conclude that be¬ 
low He there are no critical fluctuations of polarization. 
The anomalous contribution to dielectric susceptibility 
is confined to the ordered phase and is triggered by the 
spontaneous magnetic order. The situation is akin to 
conventional thermal phase transitions in improper fer- 
roelectrics, where P appears not in a spontaneous way, 
but as a result of coupling to some other parameter which 
actually undergoes criticality [1^]. Thus, the QCP in 
SUI-CU 2 CI 4 is related to the magnetic degrees of freedom 
alone, while the field-induced phase is an improper ferro¬ 
electric. 

Although ferroelectricity does not drive the QPT in 
SUI-CU 2 CI 4 , the very precise measurements of e enable 
us to extract the critical held with very high accuracy 
for this compound, and thereby elucidate the nature of 
the magnetic transition (Fig. |^. The parameter most 
relevant to the QCP is the so-called crossover exponent 
ip, which dehnes the phase boundary at T —>■ 0: 

H -HeOcT^/^. (1) 

A shrinking ht window analysis |15j of our data reveals 
that a true power-law behavior extends only up to T ^ 
300 mK, a range not probed by previous specihe heat 
studies [iniin]. However, as can be seen from the Fig.[^ 
inset, our dielectric data are dense and accurate enough 
to yield a very reliable estimate p = 0.63 ± 0.03 for that 
htting range. The anomalous result p ~ 0.34 previously 
obtained by neutron scattering [ 12 ] is likely due to the 
fact that of the hve data points obtained in this study, 
all except one did lie outside of the true power-law range. 
Our present result p ~ 0.63 is fully consistent with p = 
2/3 expected for a magnetic BEC transition [T112| . 

This finding is significant. Indeed, despite numer¬ 
ous claims to the contrary, the field-induced transition 
in most gapped quantum magnets is, strictly speak¬ 
ing, not in the BEC universality class. Instead, due to 
the magnetic anisotropy which breaks the prerequisite 
50(2) symmetry of the Hamiltonian in materials such 
as TlCuCla and IPA-CuCla, it is actually of the Ising 
type [IH] with an energy gap in the magnetically ordered 
state muD]. Even in the one known tetragonal com¬ 
pound DTN the transition is likely discontinuous 

due to magnetoelastic coupling and a spontaneous lat¬ 
tice distortion HH], showing critical exponents inconsis¬ 
tent with BEC [231 . In contrast to all these commensu- 
rately ordering materials, the 50(2) = C/(I) symmetry 
in SUI-CU 2 CI 4 is protected by its incommensurability. In¬ 
deed, the phase of the incommensurate spiral structure is 
decoupled from any magnetic anisotropy terms that are 
commensurate, and the ordered state necessarily has a 
gapless “sliding mode” |24| . 

Even though in SUI-CU 2 CI 4 the charge degrees of free¬ 
dom do not drive the QCP, the dielectric properties here 



FIG. 4. (Color online) Phase transition in SUI-CU 2 CI 4 as 
seen by nonlinear dielectric spectroscopy. Current harmon¬ 
ics In{H) (n from 2 to 6 ) induced by an applied sine voltage 
are plotted as a functions of magnetic field at T = 100 mK. 
Projections of In{H) onto Ke-H, Im-H and Re-Im planes are 
also shown. Drive voltage has amplitude F = 30 V and fre¬ 
quency u)/2'k = 551 Hz. 


are quite unusual. A careful examination of the observed 
magnetocapacitive effect reveals its amazingly nonlinear 
nature. Just varying the probing voltage leads to a dras¬ 
tic change of the amplitude of the dielectric anomaly, 
without affecting the onset of the transition or the back¬ 
ground Cq. An example of such dependence is shown in 
Fig. [^c). The peak magnitude ACmax at a fixed tem¬ 
perature monotonically decreases with a decrease of the 
voltage. Note that the electrical fields producing this 
nonlinearity are very modest [upper scale in Fig.j^c)]. 

The best indicators of nonlinear behavior are the 
higher-order harmonics in ac measurements. In our 
case of SUI-CU 2 CI 4 , we tracked the complex harmonics 
of the displacement current, induced by the applied ac 
voltage A representative data set (note the low 

T = 100 mK) is shown in Fig. Below the critical 
field the nonlinear response is zero. However, inside the 
ordered helimagnetic phase all current harmonics with 
frequencies up to 6 w show a complicated behavior in the 
complex plane as a function of applied magnetic field. At 
still higher fields, they gradually decrease and are even¬ 
tually suppressed, as is the linear magnetocapacitive ef¬ 
fect. Since the displacement charge flow is related to the 
change of sample polarization dP/dt, the appearance of 
current harmonics under a periodic voltage V sin(a;t) is 
a direct consequence of the polarization nonlinearity: 
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FIG. 5. (Color online) Nonlinear contributions to the mag- 
netocapacitive effect in SUI-CU 2 CI 4 at T = 100 mK Cn (see 
text) as a functions of magnetic field. Thick dashed red and 
dotted blue lines are real and imaginary parts of Cn and thin 
black line is the absolute value. Vertical line marks the phase 
transition field. 


P{E) — Pq + So {xi^ + X^E'^ + XaE^ + ...) . (2) 

In this expansion xi = s — 1 is the conventional first- 
order susceptibility, while x„ = are the non¬ 

linear hypersusceptibilities. The algebra required to re¬ 
store the hypersusceptibilities Xn from current harmon¬ 
ics In is summarized in Ref. |25) . To allow a quanti¬ 
tative comparison between different hypersusceptibilities 
(which are of different physical dimensions), we present 
them in the form of anharmonic capacitance contribu¬ 
tions Cn = eoXn§E^~^ [15]. The resulting Cn{H) for 
the representative T = 100 mK case are plotted in Fig.[^ 
By comparing this plot to Fig. one can immediately 
see that some of the anharmonic contributions have the 


same order of magnitude, as the linear magnetocapac¬ 
itance term. Already at a very moderate electric field 
E ~ 0.3 kV/cm, the nonlinear response is at least as 
important as the linear one. Such a giant dielectric non¬ 
linearity at a quantum phase transition of a nonelectric 
nature is an interesting phenomenon. 

A transition-related nonlinear electric response is well 
documented for proper ferroelectrics and relaxors [26] . 
but not for improper ferroelectrics. Previous studies 
of other Cu^^ based S' = 1/2 improper helimagnetic 
ferroelectrics, such as CuCl 2 [17], LiCu 202 jlHjj and 
LiCuV 04 [IHlISni, have not encountered P{E) nonlinear¬ 
ity at the transition point. A small nonlinearity was typi¬ 
cally found deep in the ordered phase as a consequence of 
a well-developed P{E) hysteresis curve. The direction of 
P in individual magnetic domains is coupled to the spi¬ 
ral chirality, which can be switched by sufficiently strong 
E, producing the hysteresis. In SUI-CU 2 CI 4 , a small his¬ 
tory dependence of magnetocapacitance above He and 
suppression of the effect by bias field [TS] indicate that 
domains also play an important role. The observed non¬ 
linearity may in principle originate from the extreme sen¬ 
sitivity of the domains to E close to the magnetic order 
breakdown. Since even tiny amounts of impurities are 
known to have a huge effect on the phase transition in 
SUI-CU 2 CI 4 [T7|, in the future it will be very interesting 
to investigate their influence on the dielectric response 
and domain mobility. 

In summary, the field-induced QCP in SUI-CU 2 CI 4 is a 
purely magnetic one, and appears to be one of the best 
realizations of magnetic BEG. The material is thus an 
“improper quantum field-induced ferroelectric.” Its un¬ 
usual dielectric response is confined to the magnetically 
ordered phase and is hugely nonlinear, even in very mod¬ 
est drive fields. 

We would like to extend our thanks to Prof. B. Bat- 
logg for generously sharing some key equipment for the 
present experiments, and to Prof. M. Mostovoy for il¬ 
luminating discussions. Special thanks go out to Dr. 
S. Gvasaliya for assistance with the instrument devel¬ 
opment. This work was supported by the Swiss National 
Science Foundation, Division 2. 

SUPPLEMENTAL MATERIAL 


I. EXPERIMENTAL DETAILS 

The sample for the magnetocapacitive measurements 
was a co-aligned mosaic of d ~ 0.5 mm-thin single crys¬ 
tals of SUI-GU 2 GI 4 was sandwiched between the copper 
plates of a capacitor, that were parallel to the crystallo¬ 
graphic (ac) plane. This assembly was placed at the sam¬ 
ple stage of Quantum Design Physical Properties Mea¬ 
surement System (PPMS) ^He-'^He Dilution Refrigerator 
insert. The magnetic field H was applied parallel to the 
a axis. It defines the spiral plane, the expected polar¬ 
ization being along b* || E. Special care was taken to 
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prevent any exposure of the sample to the atmosphere to 
avoid deterioration. Measurements were done using an 
AH2550A capacitance bridge. The capacitor plate area 
was A = 8 X 8 mm^. Due to the difficulty of precisely 
measuring the filling factor of the SUI-CU 2 CI 4 “effective 
capacitor”, we prefer to report the results as a change in 
sample capacitance C = rather than a change in 

dielectric permittivity. The values of E and e mentioned 
in the text should be considered as estimates. 

Nonlinear dielectric spectroscopy was performed on the 
very same sample assembly during the same experimen¬ 
tal run. The source voltage was applied to one of the 
capacitor plates, and SR7370 lock-in amplifier was used 
to read the displacement current harmonics. We use the 
actual value oiV = 30 V to estimate the hypersuscepti¬ 
bilities Xm but smaller voltage V = 15 V to present them 
as nonlinear capacitances C„ for consistent comparison 
with the bridge data. 

Specific heat was measured on m ~ 0.1 mg SUI-CU 2 CI 4 
sample with the help of adiabatic calorimetry option for 
PPMS ^He-^He Dilution Refrigerator insert. The con¬ 
tribution of Apiezon N grease and silver sample holder 
was measured in a separate run and subtracted. The 
measurements extend down to 200 mK; below this tem¬ 
perature the data becomes progressively plagued by huge 
contributions from nuclear spin specific heat. 

II. PHASE DIAGRAM 

To investigate the critical properties of the phase 
boundary, we have performed a windowing analysis, sim¬ 
ilar to how it has been in Ref. m- The dependen¬ 
cies of estimated critical field and crossover exponent tp 
on the cut-off temperature is shown in Fig. Around 
T ~ 300 mK the dependence of He and p on the thresh¬ 
old temperature vanishes, indicating the true power-law 
behavior. For extremely narrow fit windows the resulting 
values are slightly scattered due to the decreasing amount 
of datapoints. For a wider fitting range < 1 K, we get 
tp ~ 0.5, thus recovering the value previously found by 
specific heat measurements m- 

The observed critical filed value is different from the 
one quoted in the introduction, due to the direction of 
applied field in the present study (H||a) being different 
from that in previous neutron | 12 j and calorimetry m 
experiments (H||b). Preceding dielectric study [Hj also 
used H||a field orientation,and their phase boundary is 
apparently shifted from the values for H||b they quote 
in Fig. 4. 

III. MAGNETOCAPACITIVE EFFEGT 
A. Symmetry properties 

The spiral order in SUI-CU 2 CI 4 can be parameter¬ 
ized as (Su(i')) = Si cos(Qr) -I- S 2 sin(Qr) with Si T 


S 2 . This spiral arrangement, as it has been shown by 
Mostovoy [5] and Katsura, Nagaosa and Balatsky [4] in¬ 
duces an electric polarization P oc [[Si x S 2 ] x Q] due to 
lack of inversion symmetry between the spins and non¬ 
zero spin-orbit coupling. This effect is often called “in¬ 
verse Dzyaloshinskii-Moriya” mechanism, due to anal¬ 
ogy with the Dzyaloshinskii-Moriya antisymmetric cou¬ 
pling [3211331 • the latter case the asymmetry of the 

charge distribution between the magnetic ions allows for 
a coupling term proportional to [Si x S 2 ], which may 
stabilize a spiral phase. In the former case it is inverse: 
spiral arrangement of spins enforces spatial asymmetry 
in charge distribution. 

The distinctive feature of this mechanism is the direc¬ 
tion of polarization, lying in the spiral plane and per¬ 
pendicular to the propagation vector Q. Crystals of Sul- 
CU 2 CI 4 possess a natural cleavage plane, which is almost 
parallel to Q, and the spiral plane is fixed by the ex¬ 
ternal magnetic field, as [Si x S 2 ] || H. Hence, such a 
plane is ideally suited for the observation of magnetoca- 
pacitive effect when H T Q and E T Q, H (i. e. E 
perpendicular to the “good plane” and magnetic field 
also lies in the “good plane”). The bulk of the measure¬ 
ments was performed in this crossed fields configuration. 
However, we also adopted an alternative experimental 
setup: the same sample assembly was mounted in paral¬ 
lel fields configuration, with H || E. This means, that 
the dielectric response is probed perpendicular to the 
spiral plane. In such configuration the magnetocapac- 
itive effect is drastically suppressed (see Fig. [^: while 
at r = 500 mK a prominent peak is seen in crossed 
fields configuration, only a tiny remnant peak can be ob¬ 
served in the parallel fields configuration. This indicates 
that the non-trivial dielectric response is restricted to the 
spiral plane — exactly what should be expected for the 
“inverse Dzyaloshinskii-Moriya” mechanism. 



T 

treshold 


FIG. 6. Windowing analysis of SUI-CU2CI4 phase boundary 
obtained from dielectric measurements. 
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FIG. 7. Magnetocapacitive effect in the same sample of Sul- 
CU2CI4, measured along (E _L H) and perpendicular (E || H) 
to the spiral plane. 


B. Domain-related effects 

The magnetocapacitive effect, observed in SUI-CU 2 CI 4 , 
demonstrates a history dependence in the magnetically 
ordered phase. The example can be found in Fig. 
where the held is swept from 3 to 6 T and then back 
at rate of 10“^ T/sec. The hysteresis onset is located 
above the inhection point of the steep slope of AC{H). 
Going below the critical held “resets” the history, mak¬ 
ing the curves shown in Fig. fully reproducible over 
multiple runs. Also, note the suppression of the magne¬ 
tocapacitive anomaly by a constant (bias) electric held 
applied to the sample. All together this is a strong in¬ 
dication of domain formation above He- In zero electric 
held one may expect equal population of clockwise and 
counterclockwise spiral domains, having the opposite di¬ 
rection of electric polarization vector. Note that the spi¬ 
ral plane is hxed by external magnetic held, and hence 
the polarization of the domains must behave in Ising-like 
way. Having such domains can be a natural cause for 
the memory effect in the ordered phase. Also, it explains 
the suppression of the effect by bias electric held. Such 
held would decrease the response of P to the AC probing 
held, as it would hx the preferred direction of domain 
polarization. 
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